CLAIMS 



What is claimed is: 

1. A method for estimating at least a subset of symbol sequences from received signals 
at multiple receivers in a multiple receiver communication system wherein each of said 
multiple receivers receives a sum of signals from a plurality of transmitters sharing 
the same commimication medium, wherein said signals are formed by a sequence of 
contiguous symbols, each of which occupies a fixed time duration, and samples said 
received signals with a time interval that is the inverse of an integer multiple of said 
fixed time duration, wherein said integer multiple is no smaller than the number of said 
signals divided by the number of said receivers, wherein the channels of said signals 
have finite impulse response and are different from one another, said symbol estimation 
method comprises 

(a) sampling said received signak at multiple receivers with a time interval equal to 
a 4 fraction of said fixed time duration, wherein F, an integer is no smaller than 
T7T7-, wherein 7 is the number of said transmitters, M is the number of said 
receivers, Mt is a positive integer and generating Q sampled received signals from 
each of said received signals wherein Q is an integer no larger than Mt\ 

(b) constructing data vectors from sampled received signals by filtering them with 
channel responses between said transmitters and said receivers and multiplexing 
the outputs; 

(c) obtaining the maximum length of the channel responses of users at multiple re- 
ceivers, which is denoted as d-hl, and constructing a channel response matrix for 
each receiver, 

(d) constructing space-time correlation matrices with said channel response matrices; 

(e) constructing a second set of matrices and vectors computed from said space-time 
correlation matrices and said data vector; 

(f ) applying a recursive bisection approach on said second set of matrices and vec- 
tors to generate an output sequence, wherein said recursive bisection approach 
comprises 

forward division steps which recursively compute a first series of matrices and 
vectors initially based on said second set of matrices and vectors and yield vectors 
in said first series of vectors and matrices with dJ\Li/2\ rows where dJLi is the 
number of rows of the last element in said first series of vectors and matrices, 
until a pre-defined stopping criterion is satisfied, where [^^J means the smallest 
integer which is not smaller than x, 

backward substitution steps which first construct an intermediate solution with 
the last of the first series of matrices and vectors and recursively compute, in a 
reverse order of said forward division steps, a second series of vectors with said 
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first set of matrices and vectors until exhausting all of said first set of matrices 
and vectors; 

(g) demultiplexing the vector firom the last in said second series of vectors; whereby 
said symbols of said signals are found from said output sequence. 

The symbol estimation method in claim 1 wherein said channels are convolutions of 
physical channels and their corresponding pulse shaping filters of said channels* 

The symbol estimation method in claim 1 wherein said channels are convolutions of 
physical channels and their corresponding spreading codes of said signals in a code- 
division-multiple-access communication system. 

The symbol estimation method in claim 1 wherein said channels are convolutions of 
physical channels, their corresponding spreading codes of said signals, and their cor- 
responding pulse shaping filters in a code-division-multiple-access communication sys- 
tem. 

5. The symbol estimation method in claim 1 wherein said channel matrix comprises chan- 
nel filter responses. 

6. The symbol estimation method in claim 1 wherein the vectors in said second set of 
matrices and vectors are constructed by adding up said data vectors. 

7. The symbol estimation method in claim 1 wherein each of said channel response ma- 
trices is constructed for each said receiver by putting the channel responses between 
said transmitters and said receiver as columns in a matrix according to the user index, 
zeros are padded to said channel response matrices if necessary so that the number 
of rows in all said channel response matrices is the smallest integer multiple of of the 
maximum channel response length. 

8. The symbol estimation method in claim 1 wherein said space-time correlation matrices 
are constructed by adding up the correlation results of channel response matrices. 

9. The symbol estimation method in claim 1 wherein the matrices in said second set of 
matrices and vectors comprises three matrices Ai, Bi, and Di which are constructed 
by first finding the maximum length of channel filter responses of all users, padding 
zeros to channel filter responses whereby all of them are of d + l symbol durations 
wherein d is the least integer possible, constructing d time correlation matrices for 
each receiver with d matrices obtained from dividing the channel responses at that 
receiver into d + 1 parts, constructing d spatial-time correlation matrices for multiple 
receivers by adding up the time correlation matrices for all receivers with the same 
time shift, constructing Ai according to a biock-Toeplitz structure with spatial-time 
correlation matrices 0, • • • , d- 1, constructing Bi according to another block-Toeplitz 
structure with spatial- time correlation matrices -1, - • - , -d, and finally assigning Ai 
to Di. 



2. 
3. 

4. 
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10. The symbol estimation method in claim 1 wherein said first series of matrices and 
vectors comprises a series of tiiree dJ xdJ channel description matrices which are Ai 
the ordinary diagonal block matrix, the superdiagonal block matrix, Di the special 
diagonal block matrix, and a right-hand vector which is divided into dJxl blocks 
which are indexed from one; 

11. The symbol estimation method in claim 1 wherein said forward reduction step farther 
comprises putting said second set of matrices and vectors as the first of said first series 
of matrices vectors and computing recursively new channel description matrices from 
the latest channel description matrices and a new right-hand vector with blocks from 
linear combinations of blocks of even index and their neighboring blocks in the last 
right-hand vector. 

12. The symbol estimation method in claim 1 wherein said pre-defined criterion is that the 
length of the latest generated right-hand vector is dJ or that the superdiagonal block 
matrix in the latest of said first series of matrices and vectors becomes a zero matrix. 

13. The symbol estimation method in daim 1 wherein said backward substitution step 
further comprises solving the equations as defined by the latest of said first series 
of matrices and vectors, putting latest diannel description matrices and the solution 
as the first of said second series of matrices and vectors, assigning the blocks of the 
vector of the latest in said second series of matrices and vectors to the even blocks of 
a new vector and computing the odd blocks of the new vector from these blocks and 
corresponding right-hand vectors. 

14. The symbol estimation method in claim 1 wherein said d spatial-time correlation ma- 
trices Ti, i = 0, • • • , -rf, are constructed by 

M d+i 
m=l n=l 

where the channel response matrix at receiver m, m = 1, • • • , M, consisting of 
(d+l)J rows, where <i + 1 is said maxnnum channel length, and J coltimns, which 
are the channel responses from transmitters to the receiver m, is related to H^^ n = 
1, ••-.d + lby 



where 5{-) is the Kronecker function, v is the variance of noise for minimum-mean- 
square-error joint detection and 0 for zero-forcing joint detection. 

15. The symbol estimation method in claim 1 wherein said second set of matrices and vec- 
tors consists of channel description matrices Ai, Bx, Di and Yi where Ai is composed 
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ofdxd blocks, and the (ij) block is T-j+i if j > i, and T.j+i if j < i; Bi is composed 
of d X d blocks, and the (t, j) block is T-d+j-i if i > j, and Oyx/ if i < j; Di is AJ; 

and Ki is 53 i2< where i?i is said signal obtained from receiver i. 

16. The symbol estimation method in claim 1 wherein said second set of matrices and 
vectors are B'iArS B^A^S BiD"^ if is odd, Xi+i, and Yi, for backward substitution 
step i, t = s, 5-1, • • • , L 

17. The symbol estimation method in claim 1 wherein said forward division step comprises 
for forward division step t, t = 1, • • -,5, 



At+i 
Bi+i 
B"i+i 



-B-<A-^B< + Ai- BiA-^B'i 

-B-iA-^yi^_i +Yi^- BiA-^Yi^M, A: = 1, • • • , J 
Di+i - Ai-BiDr'B'i-B-iA-^Bi 
W^i - -B'iA-'Yij,,.2 + Yij^,.i-B\D-'Yi,L, 



if Li is odd; 



if Li is even. 



Li+i 
Yj+i^+i 



Bi - B-iA-^Bi 
= iLi 

:= -B-iA-^yi^,-i + ri,i, 



(41) 



(42) 



(43) 



18. The symbol estimation method in claim 1 wherein said pre-defined criterion is Li+i = 1 
or Bi^i = 0, 5 is set to i + 1 and the forward reduction steps stop. 

19. The symbol estimation method in claim 1 wherein intermediate solution step comprises 
solving DjXj = Ys. 

20. The symbol estimation method in claim 1 wherein said backward substitution step i, 
i = 5, • • • , 1, comprises 



Xi,2k+\ 



- A-1 



A-'Yi,i-A-'BiXi^i,i, 
A-'Yi^+i - K'B'iXi+ij, - A-^Bi.Yi+i.jk+1, 1 < 2A: + 1< Li, (44) 
D-'Yij,-D-'B'iXi+ij,,^, if Li is odd. 
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21. A method for estimating at least a subset of symbol sequences from a 
received signal at a receiver in a communication system, the method comprising: 

(a) oveisampling said received signal, thereby producing a plurality of 
received signal sample sets; 

(b) constructing data vectors from the received signal sample sets 

(c) obtaimng the maarimiTm length of the chaimel responses of users at 
multiple virtual receivers, which is denoted as d+1, and constructing a chaimel response 
matrix for each receiver; 

(d) constmcting space-time correlation matrices with said chaimel response 
matrices; 

(e) constructing a second set of matrices and vectors computed from said 
space-time correlation matrices and said data vector, 

(f) generating an output sequence by applying a recursive bisection approach 
on said second set of matrices and vectors, wherein said recursive bisection approach 
comprises: 

forward division steps M^ch recursively compute a first series of matrices and 
vectors initially based on said second set of matrices and vectors and yield vectors in said 
first series of vectors and matrices with dJ[Li/2] rows where dJLi is the number of rows 
of the last element in said first series of vectors and matrices, until a pre-defined stopping 
criterion is satisfied, where [x] means the smallest integer which is not smaller than x, 

backward substitution steps which first construct an intermediate solution with the 
last of the first series of matrices and vectors and recursively compute, in a reverse order 
of said forward division steps, a second series of vectors with said first set of matrices 
and vectors until exhausting all of said first set of matrices and vectors; 

(g) detecting symbols in the received signal by demultiplexing the vector from 
the last in said second series of vectors; whereby said symbols in the received signal are 
found from said output sequence. 



22. The method of claim 21, wherein each of said multiple receivers receives 
a sum of signals from a plurality of transmitters sharing the same communication 
medium, wherein said signals are formed by a sequence of contiguous symbols, each of 
which occupies a fixed rime dx^ation, and samples said received signals with a time 
interval that is the inverse of an integer multiple of said signals divided by the number of 
said receivers, wherein the channels of said signals have finite impulse response and are 
different firom one another. 
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23 . The method of claim 2 1 , wherein said oversampling said received signal is 
performed at multiple receivers with a time interval equal to a 1/P fraction of said fixed 
time duration, wherein P is an integer no smaller than J/MMt, wherein J is the number of 
said transmitters, M is the number of said receivers, Mt is a positive integer and 
generating Q sampled received signals from each of said received signals wherein Q is an 
integer no larger than Mt. 

24. The method of claim 2 1 , wherein said constructing includes constructing 
the d ata vectors from the sampled received signals by filtering them with channel 
responses between said transmitters and said receivers and multiplexing the outputs. 
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